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A general relation between entropy and an evolutionary superoperator is derived based on the 
theory of the real-time formulation. The formulation establishing the relation relies only on the 
framework of quantum statistical mechanics and the standard definition of the von Neumann entropy. 
Applying the theory of the imaginary- time formulation, a similar relation is obtained for the entropy 
change due to the change in reservoir temperatures. To show the usefulness of these formulas, we 
derived the expression for the entropy production induced by some dissipation in an open quantum 
system as the exemplary model system. 
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I. INTRODUCTION 

A remaining big challenge in modern physics is prob¬ 
ably to understand nniversal aspects of nonequilib¬ 
rium states in microscopic systems, and/or to develop 
nonequilibrium thermodynamics and statistical mechan¬ 
ics. Nonequilibrium phenomena are ubiquitous in nature. 
Yet, a general framework allowing their description far 
from equilibrium is lacking [l| . A defining feature of out- 
of-equilibrium systems is that they dissipate energy in the 
form of heat, leading to an irreversible increase of their 
entropy. The nonequilibrium entropy production and its 
time derivative (the entropy production rate) are there¬ 
fore two fundamental concepts of nonequilibrium thermo¬ 
dynamics [I|. The subject of time-evolution of open sys¬ 
tems with reference to the entropy production has been 
widely considered in the frame of quantum theory and 
quantum field theory, in particular, for example in the 
monumental work done by Umezawa and his school [2“ 

9 . 

Traditionally, the entropy production is used to evalu¬ 
ate the performance of thermodynamic devices: the max¬ 
imal useful work, the availability (or energy), that can be 
extracted from a given system, is reduced by the presence 
of irreversibilities, such as friction or nonstatic transfor¬ 
mations. This loss of availability is directly related to the 
mean entropy production Q . 

In recent years significant insights into the nature 
of irreversibility and entropy production have emerged, 
partly due to the need for a framework to interpret ther¬ 
modynamic processes on the nanoscale. In order to de¬ 
velop a theory for deterministic thermostatted systems, 
we take into account the dissipation function and fluctua¬ 
tion theorem [M3- Some powerful formulas such as the 
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Crooks and Jarzynski relations |lll - ll3j stand out along 
with a unifying framework for overdamped Langevin dy¬ 
namics [l3] based on a stochastic description of the first 
law of thermodynamics commonly referred to as stochas¬ 
tic energetics [13 ■ On the other hand, since the early 
ideas of Kirchoff, Rayleigh-Onsager and Prigogine 16- 
[ 13 , it was hoped, thought, and believed that nonequilib¬ 
rium regimes could be characterized by variational prin¬ 
ciples. Prigogine’s minimum entropy production princi¬ 
ple characterizes the nonequilibrium steady state as an 
extremum of the entropy production rate. For typical 
example, the stationary state in heat conduction is char¬ 
acterized by a minimum of the entropy production, com¬ 
patible with the imposed temperature distribution at the 
walls of the system [13 ■ These works indicate importance 
of obtaining non-trivial relations in analyzing nonequilib¬ 
rium systems. 


Although various available viewpoints onto the irre¬ 
versibility nature differ from each other, it is claimed that 
the central point of view is a mechanical, or possibly dy¬ 
namical, quantity can be defined throughout behavior of 
the thermodynamic entropy. Indeed, the positive aspect 
of nonequilibrium entropy production is its experimen¬ 
tal verification [ 23 , [m. Therefore, it is very important 
to find out some non-trivial relations, especially the rela¬ 
tions between entropy production and its associated evo¬ 
lutionary operator that describes the behavior of dissipa¬ 
tive systems. 

This paper is organize as follows. In Sec. im a general 
formula, which relates the entropy and an evolutionary 
superoperator for a dissipative system, is derived within 
the framework of statistical mechanics and the definition 
of von Neumann entropy in the real time formulation. 
In Sec. mil a general formula for the entropy change for 
deterministic thermostated systems is derived by intro¬ 
ducing the general theory of an imaginary-time formula¬ 
tion. In order to demonstrate a usefulness of the formulas 
obtained in this work, we consider a model for a quan- 
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turn dissipative system in Sec. lIVI Concluding remarks 
are given in the final section. 


II. REAL-TIME FORMULATION 


A general formula, which connects entropy with an 
evolutionary superoperator, is derived by using only an 
entropy operator and a causal evolution equation. In the 
microscopic picture, the state of a system is generally 
expressed by a density operator p. For every physical 
system there exists a superoperator T (the evolutionary 
superoperator), defined on the state space of the system, 
which determines the causal evolution of the state de¬ 
scriptor p via the following law: 



( 1 ) 


where the evolutionary superoperator T corresponds to 
the Liouvillian superoperator for an isolated system. The 
expression for the evolutionary superoperator T depends 
on the system considered. 

Before introducing the explicit form of T, let us first 
derive a general relation between entropy operator q and 
the evolutionary superoperator T. 

The entropy operator is introduced as a dynamical 
variable in terms of the density operator p. The entropy 
operator q is defined as 


q = —kin p , 


( 2 ) 


Comparing Eq. © with Eq. {|5j) , we obtain 


P 


pi 

0 


K(A)dA = -kTp. 


( 7 ) 


This is an exact operator identity. Taking a trace of 
Eq.(0, we can write 


K(A)dA^ =-kTrTp. 


( 8 ) 


This formula connects the entropy change rate dS/dt to 
the evolutionary superoperator T of the system under 
study. 

The equations (7) and (8) were strictly obtained from 
the causal evolution of p and the definition of the entropy 
operator q only. Noting that 


^ ^ /^\ 
dt dt \ 9t / ’ 


(9) 


we can obtain entropy S from Eq. (|S]). The last equality 
in Eq. holds since the trace operation is time inde¬ 
pendent, i.e., (9q/9t) = d(q)/dt. The similar relation 
was originally derived in the framework of thermofield 
dynamics by the present authors (23 - 1^ . 

Using Eqs. dH) and ([9]), the entropy change AS can be 
obtained in terms of the evolutionary superoperator T in 
the following form: 


AS =-k 


dtTrTp. 


( 10 ) 


where the presence of k (Boltzmann’s constant) ensures 
the von Neumann entropy denoted by S. See below. 

The expectation value of q is then given by 

(q) = Tr p q = —k Tr p In p = S , (3) 

where the statistical average of q denoted by a symbol 
(q) is the well-known “von Neumann entropy” if k is 
taken as the Boltzmann constant. We call the S defined 
in Eq. ([2|) “entropy” hereafter. From Eq. ©, the density 
operator p for the system may be expressed in terms of 
the entropy operator q: 

p = e-’i/\ (4) 


where the time dependency of the density operator p 
comes from the entropy operator q. It should be noted 
that q is indeed a function of p = p(t), i.e., q(t) = 
q[p(t)] [H. 

Taking the time derivative of Eq. (|4]), we obtain 


9t k 

where K(A) is given by 
K(A) EE 

9t 


K(A)dA, 


( 5 ) 




L9t’ k 


-H 
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L9t’ kJ 


+ 0(A3 


( 6 ) 


The relations ([5|) and ([H|) were derived by using the causal 
evolution equation ((T|) and the defining equation ([5]) for 
the entropy operator q. 

The entropy change AS induced by, e.g., the sponta¬ 
neous or irreversible processes, can be evaluated by us¬ 
ing the formula (US, which relates the entropy S to an 
evolutionary superoperator T associated to the system. 
Thus by using the obtained formula, the effect of time- 
evolution of the system due to, e.g., the spontaneous or 
irreversible process of the system can be theoretically in¬ 
vestigated by giving the explicit form of the evolutionary 
superoperator if the density operator p is given. 

Implementing Eq. (jH]) into Eq. (|3]), we receive a 
nonequilibrium distribution function at a time t in terms 
of the entropy S of the system: 


f(t) = exp(—S/k) = Cexp^— 


1 

k 



= C exp 


dtTrTp, 


( 11 ) 


where C is a normalization factor. If we adopt the initial 
condition f(to) = at t = to. Eg. (ITTl) can be 

expressed by 


f(t) 


Z V P^exp 


dt Tr Tp, 


( 12 ) 


where Z is the canonical partition function and E is the 
energy of the system at t = tg. 
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III. IMAGINARY-TIME FORMULATION 


In the imaginary time formulation (ITF), the evalua¬ 
tion of thermal averages is performed by replacing the 
time t by the imaginary time i|3, according to the formal 
analogy between the imaginary time and the inverse tem¬ 
perature first noticed by F. Bloch Time integration 
is restricted to a finite domain along the imaginary axis, 
from 0 to i(3. Thus, the time evolution in the real time 
formulation corresponds to the temperature evolution in 
the ITF. 

Let us first derive a general formula for the entropy 
production induced by the change in reservoir tempera¬ 
tures. In order to treat such a system, it is convenient 
to rewrite the formulas obtained in Sec. m by replacing 
the imaginary time it by temperature (3 (i.e., it —>■ (3). 
Here we regard the imaginary time it as temperature 
|3 (= I/lcT], and carry out the Wick rotation of the imag¬ 
inary time variable. The causal evolution equation (ED 
can then be expressed in terms of (3: 


3(3 


-iTp, 


(13) 


where H denotes the Hamiltonian of the system at the 
temperature T = l/(k(3). Equation ESD is referred to as 
the Bloch equation for the density operator of a canonical 
ensemble. 

Applying Eq. (fTOl) instead of Eq. (fT^ . Eq. (IT71) can be 
expressed and related to the entropy S by 



where H denotes the Hamiltonian. 

Hence, the generalized thermal distribution function 
f((3) is thus obtained in a similar manner as in the real¬ 
time formulation by making use of Eq. EID: 


f((3) 


exp(—S/k) = Cexp 



= C exp 


d(3{H) 


( 21 ) 


where C is a normalization factor. If we assume the ini¬ 
tial condition f((3o) = at temperature Tq = 

(k|3o)-\ Eq. m can be expressed by 


and Eq. in the Wick rotated frame is given by 

fi 


0|3 


W 

k 


Kp(A]dA, 


(14) 


where K( 3 (A) is given by 

Kp(A) (15) 

Equating Eq. E3D to Eq. EID, we obtain the exact ex¬ 
pression: 


P 


pi 

0 


K,3(A)dA = -k‘Jp. 


(16) 


In the same way to derive the formula Q, we formally 
expand K(j(A] with respect to A and then take its trace, 
and we can obtain the general formula: 



(17) 


f(|3) =Z-ie-P‘>^exp 


^ d|3(H)l , 


( 22 ) 


J Po 

where Z is the canonical partition function and E is the 
energy of the system at the initial reservoir temperature 

To- 


IV. QUANTUM DISSIPATIVE SYSTEM 

Now, we apply the formula ([S]) obtained in Sectionllllto 
an open quantum system with dissipation and consider 
how entropy is affected by the evolutionary superoperator 
associated to the open quantum system. 

A quantum mechanical transition between the states 
of the system will be considered. In the case where final 
states are continuous, we can use Weisskopf-Wigner ap¬ 
proximation . The Hamiltonian of such systems can 
be generally expressed by 


This equation relates the evolutionary superoperator T 
to the entropy change induced by the change of tempera¬ 
ture T if (3 = (kT)^^, where k is the Boltzmann constant. 
We can obtain the entropy change AS induced by the the 
temperature change in the following form: 


AS = ik 


d|3TrTp. 


(18) 


Now we consider a canonical ensemble. In thermal 
equilibrium, the density operator is given by p ~ 

This unnormalized density operator satisfies the equation 


3(3 


= -Hp, 


(19) 


H = Ho+ir, (23) 

where Hq and E are Hermitian operators, denoting the 
Hamiltonian of an open system and the system-reservoir 
interaction, respectively. 

In this model, the time-evolution of the density oper¬ 
ator follows the Liouville-von Neumann equation: 

= “^[Ho, p] + ^{E, p} = Tp , (24) 

where the curly brackets indicate the anticommutator de¬ 
fined for any operators A and B: {A, B} = AB -|- BA. If 
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p commutes with Hq (i.e., [p, Hq] = 0), the entropy pro¬ 
duction rate can then be obtained from the formula (jHl) 
as in the form: 


dS /\ 
dt ^ W/ 



(25) 


where P is an operator representing dissipation (damp¬ 
ing). We can evaluate the entropy production rate by 
giving a specific form (microscopic description) of P and 
thus obtain the P-related information theoretically. The 
obtained formula enables us to obtain some information 
about the entropy directly from the evolutionary super¬ 
operator “J in the equation of motion for the density op¬ 
erator p. 


V. CONCLUDING REMARKS 

A microscopic theory of entropy was developed and a 
general formula connecting entropy to the evolutionary 
superoperator associated with a system was derived. We 
also derived some general relations by making use of only 
the causal evolution of the density operator and the def¬ 
inition of the von Neumann entropy in the framework of 
quantum statistical mechanics. Making use of the rela¬ 
tions, it is easy to obtain the entropy change rate directly 
from the evolutionary superoperator T. Also we showed 
that by making use of these relations, the nonequilib¬ 
rium distribution function for the system can be derived 
by taking the expectation value of the entropy produc¬ 
tion, see Eq. m in the real time formulation. It is worth 
stressing that that the distribution function has the form 
which is not an expectation value of exponential func¬ 
tion but an exponential function of the expectation value. 
Consequently, we have obtained a specific form of the dis¬ 
tribution function of the concerned system. 

In the imaginary-time formulation we obtained the re¬ 
lation (20) from the Bloch equation and the statistical 


entropy in a similar manner to the real-time formulation 
discussed in Sec. 2. We also derived the generalized ther¬ 
mal equilibrium distribution function from Eq. (1^ in the 
imaginary time formulation. The distribution function 
in imaginary-time formulation recover the conventional 
thermal equilibrium distribution function by taking the 
lowest-order approximation for the entropy production. 

The relations obtained in Sec. 2 can be applied to an 
open quantum system. In order to show the applicability 
of obtained formulas, we considered the case where the 
interaction causes dissipation. We obtained the formula 
for the entropy production due to the dissipation arising 
from a specific interaction by giving the specific interac¬ 
tion Hamiltonian. In other words, if we find the micro¬ 
scopic description of the interaction Hamiltonian P, all 
we need to do is to utilize the relations (jH]) and (HOD, and 
we will have the entropy production from which we can 
obtain some information about the interaction through 
the entropy of the system for the open quantum system. 
In the standard approach, entropy production is achieved 
through nonequilibrium thermodynamics. In our theory, 
however, we can directly evaluate entropy production by 
employing the formulas given in this paper. 

Finally it is worth mentioning that a possible extension 
of our work would be a construction of the microscopic 
foundation of nonequilibrium thermodynamics. In par¬ 
ticular our results would have the possibility to obtain 
microscopic representations of the system-reservoir inter¬ 
actions i.e. heat. Our new formulation presented in this 
article enables us to calculate the entropy production di¬ 
rectly from system-reservoir (environment) interactions. 
This means that it allows us to conjecture what type 
of interactions play an essential role for frictional effects 
on the system through application of a conjectured form 
(microscopic description) of the interaction Hamiltonian 
r. It goes beyond the standard derivation of the second 
law of thermodynamics. 
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